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Abstract 

A theoretical description of the dissociative recombination process for the HCO + ion suggests 
that the nonadiabatic Renner- Teller coupling between electronic and vibrational degrees of freedom 
plays an important role. This finding is consistent with a recent study of this process for another 
closed-shell molecule, the Hg~ ion, where Jahn- Teller coupling was shown to generate a relatively 
high rate. The cross section obtained here for the dissociative recombination of HCO + exhibits 
encouraging agreement with a merged-beam experiment. 

PACS numbers: 34.80.Ht 34.80.Kw 34.80.Lx 
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I. INTRODUCTION 



The dissociative recombination (DR) of small molecular ions that collide with electrons 
plays an important role in interstellar diffuse and dense clouds. It is well known that these 
clouds constitute building material for new stars. The importance of DR cross sections 
as parameters in star formation models, for such astrophysically relevant ions as H^" and 
1C0 + , is one of the reasons why DR has been extensively studied in laboratory experiments 
l|. For diatomic molecules the process is well understood and described theoretically 
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Until recently 
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12j though, theory was unable to model the DR 



process in triatomic ions, except in cases where a neutral dissociative state crossed the ionic 
Born-Oppenheimer surface in the Franck-Condon region, generating a rapid rate |9']. In 
triatomics, a key complication is the fact that multiple vibrational and rotational degrees of 
freedom must generally be taken into account. And in addition to the greater computational 
burden of treating more dimensions quantum mechanically, the addition of these new degrees 
of freedom can also lead to new conceptual issues, related to the degeneracy of vibrational 
or rovibrational levels in certain triatomic ions (e.g., in both and HCO + ). This can 
cause an intrinsic instability of the corresponding neutral molecules that causes them 
to distort away from the symmetric geometry. 

Consider an incident electron that interacts with a closed-shell triatomic ion having a 
degenerate vibrational mode. If the symmetry group T of the resulting neutral complex 
has at least one degenerate irreducible representation, the electronic partial wave compo- 
nents with angular momentum I > typically contain at least one such representation, 
whereby the corresponding electronic states of the neutral molecule are degenerate in the 
clamped-nuclei approximation. Due to the Jahn- Teller theorem Q, Q] , first formulated 
by Landau [13I . Ilq . this degenerate electronic state strongly interacts with a degenerate 
vibrational state of the T group. The interaction leads to either a quadratic (Renner- Teller) 
[nl llfl or linear conical (Jahn- Teller) intersection, where "quadratic" or "linear" refers 
to the dependence of the adiabatic eigenvalues on the symmetry distortion coordinate; e.g., 
the electronically-degenerate linear molecule being considered in the present article displays 
a quadratic behavior near the symmetry point, as a function of the bending angle 9. In 
electron-induced reactions, electron capture can be followed by dissociation or autoionization 
of the recombined system. (In some rare cases, recombination could also occur radiatively 
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following electron capture, but this is rarely important in molecular systems having an open 
dissociation channel and it will be ignored throughout this study) Since the degenerate elec- 
tronic state causes the instability of the symmetric configuration, the recombined molecule 
quickly distorts away from the symmetry point to remove the degeneracy, after which the 
autoionization channel typically becomes energetically closed, and the molecule eventually 
finds a pathway along which it can dissociate (provided there is an electronic state of the 
molecule that is open for dissociation). In a system for which no direct crossing of a neutral 
dissociative state with the ground state of the ion is present, this indirect mechanism can 
become dominant, and it can still produce large DR cross sections. This is expected to be 
particularly true for molecules containing hydrogen. 

Work by some members of t he p resent collaboration have recently found that this is what 
happens in Hjj" + e~ collisions 12] . For HCO + , a closed-shell linear ion in its ground state, 
the picture appears to be similar. The lowest doubly-degenerate vibrational E states are 
coupled to the electronic states Eiinpir) of the neutral system through the Renner- Teller 
interaction, resulting in a large probability for recombination. 

Although Renner- Teller coupling is a well-known phenomenon in spectroscopy, the ef- 
fect of Renner- Teller coupling on electron-molecule collisions has not been studied until 
recently. McCurdy et al. have made an important contribution to understanding the 
role of Renner- Teller physics in collisions between an electron and the CO2 molecule. In a 
completely ab initio treatment, the authors employed a time-dependent framework to de- 
scribe the nuclear motion of the negative ion CO^~ formed during the collision. The nuclei 
move on the two-component 2 A\ and 2 B\ electronic complex potentials corresponding to the 
doubly-degenerate 2 I1 M molecular state at linear configuration. The two components interact 
through the Renner- Teller coupling mechanism. 

The aim of the present study is to understand the mechanism of DR in HCO + and, 
in particular, to interpret its large measured rate. This ion was first detected in space 
and then synthesized in the laboratory j^J. Experimental measurements of DR in HCO + 
were realized in a number of merged-beam, afterglow-plasma, and storage ring experiments 
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In Ref. |31( we have presented the potential energy 
surfaces for the ground molecular state of HCO + and the lowest states of HCO up to principal 
quantum number n — 4, calculated as functions of all three vibrational coordinates. The 
potentials were then used to estimate the DR rate. Since the splitting between the two 
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Eiinpir) electronic states for the principal quantum number n = 3 was found to be small, 
it was suggested that the Renner- Teller non-Born-Oppenheimer coupling should not play a 
significant role in DR of HCO + , therefore, it was not included in the model of Ref. |31|. 
However, since even a small coupling between degenerate states can in principle result in a 
large DR probability, we investigate here the role of Renner- Teller coupling in HCO in more 
detail, and find that we must revise that initial expectation. 

The article is organized in the following way. In section |H] we construct the model 
Hamiltonian from the results of the ah initio calculation. The Hamiltonian includes the 
nonadiabatic Renner- Teller coupling and depends parametrically on the nuclear configura- 
tions. Section ITTT1 describes the quantum defect approach to obtain the reaction matrix from 
the Hamiltonian. It also describes how the DR cross section is obtained from the reaction 
matrix. Section[lV]presents and discusses the results of the calculation: the DR cross section 
and the thermally averaged DR rate. Section IS considers autoionization in detail. Autoion- 
ization is a process that competes with DR and it decreases the DR cross section. Section 
IVTI compares the present results with a previous theoretical study j^]] • Finally, section IVHI 
discusses our conclusions. 



II. MODEL HAMILTONIAN INCLUDING THE RENNER-TELLER COUPLING 

The starting point of the present theoretical approach is the molecular Hamiltonian H of 
HCO, which we represent as H = H ion + T e + H int , where H ion is the ionic Hamiltonian, T e is 
the kinetic energy of the incident electron, and H int describes the electron-ion interactions. 
We assume that the ion is in its ground electronic state. H ion and H int depend on the four 
internuclear coordinates Q = {Rch, Rco-i^-i^p}, where Rch and Rco are the C-H and C-0 
internuclear distances, 9 is the bending angle, which is zero for linear configurations. The 
electronic energies are independent of the angle if, which represents the azimuthal orientation 
of the bending. 

We assume that the incident electron is initially captured in one of the neutral Rydberg 
states, which can be approximately characterized by the principal quantum number n, the 
orbital angular momentum I, and its projection on the molecular axis A. In our model we 
have included only the npir~ l ,npir +1 , npa, nsa and nda states. The symbols ±1 imply two 
opposite- sense electronic angular momenta associated with the different signs of A. Linear 
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combinations of the nn^ 1 states give states npir' and npir", symmetric and antisymmetric 
with respect to reflection in the plane containing the molecular axis. The other three states 
are symmetric. These states were established in Ref. 3jJ| to be the most important for 
DR, since they exhibit the largest dependence on the vibrational coordinates. Therefore, 
H int +T e assumes a block-diagonal form with an infinite number of 5 x 5 blocks corresponding 
to n = 2, 3, • • • , oo, in addition to the continuum. Since the only off-diagonal couplings 
included in our analysis are those among the three np states, for simplicity we specify only 
these states in the formulas below. 

The np-block of H int + T e in the basis of the npn' 1 , npa and npir +1 states has the form 
3J 

E v Se* 7 e 2i< ^ 
H int {Q) = 5e-^ E a Se* , (1) 

where E a and E n are the electronic energies of the npa and np7r ±l states at the linear ionic 
configuration; 5 and 7 are the real, non-Born- Oppenheimer coupling elements. We denote 
both the (me — nir) and (mr — no) couplings as Renner- Teller (RT) couplings, whereas in 
some previous studies, only the former is denoted by this term. The couplings S and 7 
depend on Rch, Rco, an d 9. The diagonalization of [T e + H int \(Q) is accomplished by the 
unitary transformation matrix U 
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with the abbreviations 



w± = >/(l±A/w)/2; A = {E a - E n - 7 )/2; 



w 



VA 2 + 25 2 . 
When diagonalized, the Hamiltonian becomes 

U\H int + T e ]U = diag{V^, V,, V^} 

where 



(3) 



(4) 



K» = K-7, V a = (E a + E n + 1 )/2 + w , 
V n > = (E a + E v + -y)/2-w. 



(5) 



The adiabatic potential energy surfaces V n iy^(Q) are known from ab initio calculations 
(see Ref. for a detailed description). From Eq. (j3J) we obtain 

7 = E n -V^, w = (V a - VV)/2 • (6) 

Therefore, the matrices U in Eq. (j2J) and H int in Eq. JTJ) are obtained from K-/ )7 r",o-(Q)- 



III. REACTION MATRIX AND CROSS SECTION FOR DISSOCIATIVE RE- 
COMBINATION 

As in our previous DR studies ^j, U, 3|, we employ multichannel quantum defect th eory 
and need to construct the reaction matrix K, which is related to the potential as in Ref. 1331. 

is not the i^-matrix of MQDT, 



which includes an additional factor of — 7T.) First we introduce the diagonal quantum defect 
matrix /i, whose nonvanishing elements /i, are = — l/[2(n — /ij) 2 ], where i is the electronic 
state index. Therefore, the diagonal form of the reaction matrix is directly obtained from 
the ab initio calculation. 

The coupling 5 in Eq. Qis linear with 9 for small 9. Correspondingly, the splitting between 
adiabatic electronic energies V n n and V n > is quadratic with respect to 9, and, therefore, the 
splitting between quantum defects /v and is also quadratic. This can be seen in Fig. 
Q where /!„-» is almost constant because it has a different symmetry and is unchanged due 
to the Renner- Teller effect. The quantum defect /v is approximately quadratic due to the 
Renner- Teller effect until 9 m 50°, beyond which the linear approximation for 5(9) is not 
adequate. In a preliminary calculation we used such a linear dependence of 5(9), because we 
have found that we can use formulas Eqs. El El El El and El and fit the model Hamiltonian 
accurately to the ab initio energies without the additional approximation of linearity of 5(9). 
The result of the preliminary calculation is consistent with the present calculation using the 
accurate fit to the ab initio energies. 

The next step is to bring the reaction matrix into the representation of the Hamiltonian 
of Eq. ^ Since the transformation matrix U is known from Eq. EJ the K matrix in that 
representation reads 

K = U tan(7r/i) . (7) 
The matrices /t, H int + T e , K, U are diagonal for the nscr and ndcr states. We used quantum 
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FIG. 1: (Color online) The quantum defects of the electronic states included in the present study 
as a function of the bending angle for fixed Rch = 2.0 a.u. and Rco = 2.0877 a.u. Labeling 
the states with electronic momentum projection A is an approximation and for large values of 9 
it is not appropriate. However, the four sa, pa, da, and pn' states are always uncoupled from 
pir" states for nonlinear geometries. In the absence of the non-diagonal matrix elements in Eq. ^ 
the pn' and pn" states would have the same quantum defects. The Renner- Teller parameters are 
completely determined by the geometry dependence of pn and pa defects (Eq. EJ). 

defects from Ref. with n = 4 for nsa and n = 3 for the other states. Figure gives the 
quantum defects as functions of 9, with Rch and Rco fixed at the ionic equilibrium values. 
Once the reaction matrix Ka' is obtained, the DR treatment is along the lines of Refs. 
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34j with the following differences. We choose here the Rch distance as the dissociative 
adiabatic coordinate. This assumption is justified because experimentally the H+CO chan- 
nel is largely dominant at low energies. It has been long known that the CH+O channel is 
endothermic by 0.17 eV and could have an energy barrier, while the OH product is not ob- 
served (see Ref. j2l| and references therein). Recently, branching ratios have been measured 
3o| for the DCO + dissociative recombination, confirming that the D+CO channel is by far 
the dominant, with a ratio of 0.88. We keep the Rco distance frozen at its ionic equilibrium 
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FIG. 2: (Color online) Several adiabatic potential curves for the HCO + ion as a function of the 
Rch distance for different projections m v of the vibrational angular momentum. The curves of the 
same color correspond to the vibrational states with the same quantum number but different 
L. The energies of the lowest vibrational states with = and m v = ±1 are shown by two 
horizontal lines. 

value. This assumption was made to simplify the treatment but it could result in an un- 
derestimation of the DR cross section. In order to account for the CO vibration one might 



use the hyper-spherical coordinates similar to Ref. 31], or else normal mode coordinates 
as in McCurdy et al. [19]. Although the CO bond is frozen, the Renner- Teller coupling 
physics is included in our model; we show below that it does increase the DR cross section 
significantly. We will briefly address below the possible influence of the CO vibration. 

Consequently, for every Rch distance we calculate matrix elements of the "large" reaction 
matrix JC(Rch) 

KjARch) = {$ mv ,L\K i4 ,(Q)\$m>,,v)0,<p, (8) 

where l{Rch\ @, V 9 ) are vibrational wavefunctions of HCO + parametrically dependent 
on Rch, while Rco is frozen. They are eigenfunctions of H ion with eigenvalues l (Rch)- 
Each index j in fC can be represented as {i,m v ,L}, where m v specifies the projection of 
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TABLE I: The table demonstrates the accuracy of the vibrational energies obtained in the present 
study neglecting the CO vibration with the exact calculation from Ref. 3fij |. Since the CO bond 
is frozen, we don't provide energies for excited v% modes. The overall error is about 12 %, which 
translates into about 25 % for vibrational wave functions. The energies are given in cm -1 . 



the vibrational angular momentum on the molecular axis, while the index L distinguishes 
different vibrational states with the same m r The rotation of the whole molecule is not 
considered. Fig. El shows the vibrational eigenenergies L as a function of Rch- 

It is informative to compare the vibrational energies obtained in our simplified approach 
with the exact calculation accounting for all four vibrational coordinates puzzarini96. We 
cannot compare the absolute values of the energies because the CO stretch is frozen, however, 
we can compare the energies of a few vibrational levels {viV l 2 , v^} with respect to the energy 
E = of the ground vibrational level {00°0}. Table H] compares the vibrational energies 
obtained in this study with accurate results from Ref. [36]. The overall accuracy of the 
present calculated vibrational levels is about 12 % or less. Since the vibrational wave function 



error is of order \/0.12 « 0.35, a conservative estimate of the error in our calculations of the 
final cross section would be of order 70%. 

The matrix JC(Rch) is used to obtain the potential curves U a (Rcn) of HCO as described 
in Ref. Q]. The curves in general have non-zero autoionization widths T a (RcH)- 

The potential curves U a (Rcn) and their autoionization widths T a (R C n) are then used to 
calculate the DR cross section. The cross section is calculated in a manner similar to the 
procedure described in Refs. [ill 3], with some adaptations. Specifically, we start from 
Eq. (5.19) of Ref. for the cross section of the process of dissociative attachment of 
the electron to a neutral molecule, which applies equally to the dissociative recombination 
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process considered here. The resulting approximation adapted to this situation reads (see 
Eq.6 of Ref.H): 

2n 2 T a (R C H) I +/ D m2 -p(E) (n\ 

aDA = l^WiR^)\ Wo{RcH) ^ (9) 

Here e _p ^ < 1 describes the survival probability, which could be less than unity if the 
system has a substantial probability to autoionize before it dissociates; k Q is the asymptotic 
wave number of the electron incident on vibrational level o of the target molecule. The 
distance Rch is a dissociation coordinate understood to be evaluated at the Condon point, 
which depends on the total energy as well as the initial target vibrational level o; xt {Rch) 
is the initial vibrational wave function of the ion. The above formula is appropriate for 
the capture into a resonant potential curve that is energetically open for direct adiabatic 
dissociation. If the corresponding resonant state is closed, i.e. bound with respect to dis- 
sociation, it requires modification. To this end, we adapt Eq. (4.2) of Bardsley j^l to our 
present situation involving indirect DR. The indirect process proceeds via capture into a 
bound (typically Rydberg) state, which eventually predissociates. For a case involving a 
single incident electron partial wave (e.g. pir', for definiteness) and a single ionic target 
state, the fixed-nuclei autoionization width of the resonance potential curve will be denoted 
Fo,(Rch)- Once the vibrational motion in the resonance potential is quantized into a vibra- 
tional resonance level, with radial wavefunctions x ms {Rcii), it should be remembered that 
only a subset o of the target vibrational levels v will be energetically open, at any given 
total energy E. Each resulting quantized resonance acquires a partial autoionization width 
T a o i into an open vibrational channel o', which is given approximately by 

poo 

Ta,o> = \ xt>(RcH)Vr a (RcH)X IeS (RcH)dR C H\ 2 ■ (10) 

Jo 

The sum of these is then the total resonance autoionization width (within this approxima- 
tion, neglecting Rydberg level perturbations of the MQDT type, which could sometimes 
produce complex, non-isolated resonances), 



a, tot 



while the total linewidth T of this quantized resonance also includes its predissociation 
partial linewidth, T^: r = + r a tot . In this notation, the contribution to the total DR 
cross section will be the following, if only the ground vibrational level o of the target is 
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populated: 

rcs/7-iN ^7r 1 Tg :0 Td 

" (£) = lf^(B-B ra p + ip' (12) 

Our model does not account for vibrational motion along the CO bond and, therefore, 
positions of resonances cannot be compared directly with the experiment. Moreover, the res- 
olution in the available experiments is usually not sufficient to resolve individual resonances. 
This suggests that we should average the cross section over energy. Each energetically-closed 
ionic bound vibrational state c generates a Rydberg series of resonances e„ c , numbered with 
the effective principal quantum number n c . The cross section averaged over the energy 
interval between two resonances with energies e n and e n+ i is given by 

(-(E)) = 1 / a(E>)dE>, A n+1 = ^-^. (13) 

Extending the limits of the integral to infinity, one can easily obtain the following: 

where the A = A n + A n+i 1/nj?. Note that the quantity in parentheses here approaches 
a constant at sufficiently high effective principal quantum numbers n c in the relevant closed 
channel c, because each partial and total width in the parentheses of this formula should 
become proportional to A in this limit. The total cross section is then calculated by summing 
up the average contribution from all Rydberg states (a sum over closed ionic channels). In 
the limit where T d » r a tot , this gives 

^ = T?E \(7r(RcH)\VT^H)\xt(RcH))\ 2 nl (15) 

The projection M = + A of the total angular momentum on the molecular axis O-C 
is a conserved quantity in our model. We calculate all the resonances and the cross section 
for all allowed values of M. Since we are considering only sa,pcr, da,pir', n" quantum defects 
and the initial ion is in its ground state = 0, the only possible values of M for the final 
state are and ±1. Resonance energies and widths are of course dependent only on \M\. 
Finally, we obtain for the cross section: 

(a total ) = (a M=0 ) + 2(a M=1 ) . (16) 
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FIG. 3: (Color online) The figure shows the present theoretical DR cross section (dashed and dot- 
dashed lines) for HCO + as a function of the incident energy E. The initial vibrational state is the 
ionic ground state for the dot-dashed curve and the first excited state for the dashed curve. The 
^ Q < OT Symb o ls , md pr ev ioUS t — Q (tun S o, id Me, cross S ec t , nS a re 
also shown for comparison. The theoretical curves include the averaging over the electron energy 



distribution according to the procedure described by Eq. (2) of Ref. 
meV. 



3J with AE ± = AE\\ = 3 



IV. RESULTS AND DISCUSSION 



We have calculated DR cross sections for two different initial vibrational states xt(Rcif) 
of the ion: the ground state v = (which can be viewed as {00°0}) and the first excited 
state v = 1 (which can be viewed as {01 1 0}) whose energies are represented in Fig. El with 
horizontal dotted lines. The energy difference between the two states is about 0.1 eV. Figure 
El shows the two cross sections with dot-dashed (v = 0) and dashed (v = 1) lines. Fig. El 
shows that the DR cross section depends strongly on the initial vibrational state. A similar 
strong dependence on the initial vibrational state has also been observed in experiments 
with other molecular ions. Assuming that the the initial experimental vibrational state 
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distribution is in thermodynamic equilibrium, it is necessary to average the cross section 
over the thermal distribution to compare with the experiment. The first excited vibrational 
state has vibrational symmetry different from the ground state. The vibrational angular 
momentum for the states are m v = 1 and 0. Thus, the deexcitation process may be too slow 
to reach equilibrium, resulting in a vibrational temperature that is higher than the electron 
temperature (300K in the experiment of Ref. j^lj ). The vibrational temperature in the 
experiment of Le Padellec et al. j^lj is not known. If it was in equilibrium with the electron 
temperature at 300K, the contribution from the excited vibrational states would be small. 
But for a larger vibrational temperature, for example 1000 K, the averaged cross section 
would be about a factor of 1.5 larger than the cross section for the ground vibrational state 
in Fig. El The fact that in several different experiments the measured DR rate ranges over 



values from 0.65 - 3 x 10" 7 cm 3 /s at 300 K (see Ref. |22j, |23j, |2J, [M M, UA M, [M \M and 
Fig. HJbelow) might conceivably derive from differences in the initial vibrational populations. 
Determination of the actual experimental vibrational distribution and/or controlling it in 
HCO + could be an important step in understanding DR in small polyatomic ions. 

The present theoretical DR cross section is approximately a factor of 2 smaller than 
the experimental data. Assuming that this reflects a limitation of the present theoretical 
description this might derive from our approximation that freezes the CO bond length. 
However, the expected error of order 70% caused by our adiabatic approximation in the Rch 
coordinate should also be kept in mind when assessing the implications of this discrepancy. 
Still, if the CO bond is allowed to vibrate, which makes the ion more floppy, the probabilities 
of capturing and predissociation will presumably be increased. Quantitatively, releasing 
the CO bond increases the density of HCO resonance states that can be populated in the 
electron-ion collision, and correspondingly, the sum in Eq. E3is expected to increase. 

In afterglow plasma experiments, the measured observable is the DR rate thermally 
averaged over the kinetic energy distribution of colliding electrons and ions. This thermally 
averaged DR rate a(kT) is obtained from the DR cross section shown in Fig. El according to 
Eq. (7) of Ref. The resulting theoretical rates (thermally averaged over the electron 

energies, for each of the two initial vibrational states) are shown in Fig. (J3j) and compared 
with available experimental measurements. 

Although the theoretical DR cross section is still smaller by a factor 2 than the exper- 
imental results of Le Padellec et al. j^j], we find that inclusion of Renner- Teller coupling 
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FIG. 4: (Color online) Theoretical and experimental DR thermal rates for HCO + . The dot-dashed 
and dashed lines are the theoretical rates obtained for the ground and first excited vibrational states 
of the ion. The experimental rate (solid line) is obtained from the experimental cross section, Ref. 
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The results from a number of other experiments with HCO + 
are also shown. 
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significantly increases the theoretical DR cross section. The previous theoretical study by 
Larson et al. [3^ that omitted Renner- Teller coupling, but which did account for the CO 
bond vibrations, gave a DR cross section significantly smaller than the present result. This 
is evidence for the important role of Renner- Teller coupling in the HCO + DR process. 

Since the Born-Oppenheimer potential surfaces are the same for HCO (HCO + ) and DCO 
(DCO + ) we did a similar calculation for DR in DCO + . The resulting DR cross sections for 
the ground and first excited vibrational level are very similar to the ones presented in Fig. 
0] but smaller by a factor of 1.5. A similar dependence of the DR rate on the isotopologue 
masses was found in theory and experiment for the ITt, H2D + , D2H + , and Dj ions I2I Isoj] - 
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V. ESTIMATION OF THE EFFECT OF AUTOIONIZATION AND CO VIBRA- 
TION ON THE CROSS SECTION 



We should mention some approximations made in the present study that can affect the 
theoretical DR cross section. First, our treatment as implemented here has not accounted 
for the possible autoionization after the electron is captured by the ion, before the neutral 
molecule has time to predissociate: in Eq. (0), the survival probability e~ p ^ was set to be 
1. The inclusion of autoionization decreases the calculated DR cross section. The effect of 
autoionization on the DR cross section can be estimated as following. Figure El shows the 
calculated resonance curves U a (Rcn) as described above and in Ref. [10(. Above the lowest 
ionic curve each of these potential curves has in general a non-zero width T a (RcH), which 
represents the adiabatic autoionization rate. If there is more than one open ionic channel 
\j) for a given U a (RcH), the neutral molecule can decay into the i-th channel with a partial 
width Ya\RcH)i where Ya\Rcn) = ^a(RcH)- Once the electron is captured, the system 
can evolve into two competing pathways: resulting in either autoionization or dissociation. 
The relative probabilities per unit of time for the two processes provide an estimate for the 
survival probability e~ p<yE \ 

The largest partial autoionization widths are in general the widths of the resonance curves 
Uc{Rch) belonging to Rydberg series associated with a few nearest closed ionic channels. 
Such resonance curves are very similar in shape to their parent ion potential curves, when 
the principal quantum number n is high. Consequently, in the energy range of interest, these 
curves U a {RcH) are closed with respect to adiabatic dissociation. The dissociation can occur 
only through coupling to true dissociative states (i.e., the process denoted predissociation 
in molecular spectroscopy). The predissociation probability can be estimated using the 
Landau-Zener model. When the neutral molecule is vibrating along the U a (RcH) curve, 
every time it passes through an avoided crossing with another curve U a '(Rcn), it can jump 
to the corresponding state \a') via an adiabatic transition. We view such an adiabatic 
transition as the pathway to predissociation. The probability for adiabatic passage through 
an avoided crossing is given by the Landau-Zener formula [13] 

( 27r « 2 / ™ \ 27ra 2 / m 
Pa ' a = 1 - 6XP \ ~AF~ \ 2AEj W AF V 2AE (17) 

where a is the non-diagonal coupling element between diabatic states, i.e. the states that 
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2.1 

R rT r (a.u.) 



2.2 



FIG. 5: (Color online) Resonance potential curves U a {RcH) of the neutral molecule states having 
magnetic quantum number M = 1. The potential curves were calculated in the energy range where 
at least one electron-ion channel U^~(Rch) (for the given Rqh) is open. Thus, the lower bound of 
the calculated curves corresponds to the ground potential curve Uq (Rch) of the ion. The potential 
curves of the ion are also shown by dashed thick lines (see also Fig. |2J). The Rydberg series are 
crowded just below the ionic potential curves. For clarity, visible gaps between the Rydberg series 
and the ionic curves have been intentionally introduced artificially, to show the behavior of the 
resonance curves belonging to higher ionic channels. The horizontal line represents the energy of 
the ground vibrational level. 

would cross if the coupling was absent. Numerically, a is equal to half of the adiabatic 
potential curve splitting at the avoided crossing. AF is the absolute difference in slopes 
(net classical force) between crossing diabatic potential curves, AE is the classical kinetic 
energy at the avoided crossing, m is the reduced mass of the system. The resonance curves 
shown in Fig. demonstrate numerous avoided crossings. 

Table |H] conveys an idea for the order of magnitude of the parameters in Eq. El This 
table shows two examples of avoided crossings that lie well below the energy at which the 
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Rch a.u. 


Energy (eV) 


a (a.u.) 


AF (a.u.) 


AB(a.u.) 


Pa', a 


2.056 


0.112 


1.18 x 10~ 5 


0.0041 


0.0059 


0.9 x 10" 4 


2.11 


0.144 


2.9 x 10" 5 


0.0111 


0.0047 


2 x 10~ 4 


1.93 


0.34 


4.0 x 10" 5 


0.0036 


l-37xl0" 4 


14 - 83 x 10~ 4 


2.02 


0.283 


1.0 x 10" 5 


0.0022 


l-37xl0" 4 


1.2-7.2xl0~ 4 



TABLE II: The table demonstrates Landau-Zener parameters and probabilities for typical avoided 
crossings. The first two examples correspond to avoided crossings situated deeply below the energy 
of the ground vibrational level. The third and fourth examples corresponds to energies around and 
above the ground vibrational level. Thus, AE could be very different for the last two examples 
giving different probabilities. In general, for any total energy of the system there are always some 
avoided crossings with small AE at the corresponding turning points. 

electron can be captured. At such large values of AE, the Landau-Zener probability 10~ 4 , 
depends only weakly on the electron energy. The two other examples are taken from the 
region around the left turning point where the velocity of motion could be small or large 
depending on the energy of the electron. In cases like this one, we calculated the probability 
P a i : a for a range of values of the kinetic energy AE. There are many such avoided crossings 
that occur near left turning points and comparatively fewer at right turning points. 

During half (t!/ 2 ) of an oscillation period, the system goes through n c avoided crossings. 
We estimate that n c is about 10 for a typical resonance curve. Therefore, the total predisso- 
ciation probability Pd is somewhere in the range 0.001-0.01 . In the above estimation of the 
predissociation probability we have assumed that only two states interact at each avoided 
crossing. However, in reality many of the avoided crossing cannot be considered as strictly 
two-state ones: three or more states \a) can participate. The probability of multistate 
crossing might conceivably give a larger predissociation probability. 

The autoionization probability Pi on during time T1/2 is P.i on ~ T1/2L, where V is the 
total autoionization width. The largest Y is about 10~ 5 a.u., ti/ 2 ~ n/u, where ui/2 is 
the frequency of oscillations in the ground vibrational level, hu/2 = 0.007 a.u. Thus, 
Pion ~ ttL/c<j rs 0.002. This is presumably an upper bound on the autoionization probability. 

As one can see, the autoionization probability could be competitive with predissociation 
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0.13 



7 



0.116 



0.118 - / 



1.9 



1.92 



1.94 




R rH (a.u.) 



FIG. 6: (Color online) A part of the spectrum of Fig. [^showing the interval of Rch corresponding 
to the left turning point of the vibrational motion. For convenience, the energy of the ground 
ionic curve Uq (Rch) was subtracted from all the data. In the figure, curves with positive and 
small negative slopes correspond to Rydberg states converging to the nearest ionic thresholds. The 
probability to capture the electron tends to be high for such states, diminishing however at high 
principal quantum numbers as n~ 3 . 

in our model if we take the smallest ratio Pd/Pion ~ 0.5, but it is more likely that the 
ratio is of the order of 2-3 and therefore, more favorable to predissociation. In the previous 
theoretical study of DR in Hg" n|, the ratio was significantly larger. An important difference 
with the Hjj" study is the two degrees of freedom taken into account when adiabatic potentials 
were calculated. As a result, the density of resonance curves was much higher (factor 
20-40). The main contribution to the high density of resonance curves was from the states 
with small principal quantum number corresponding to highly excited ionic channels. This 
increased significantly the probability of predissociation. In fact, in the case of HCO, the 
vibrating CO bond should also increase the density of steep resonance curves in Figs. 03 
and El This is because if the vibration along CO is quantized, each of the curves shown in 
Fig. El will produce a series additional ionic curves with different CO quanta. Therefore, 
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the number of ionic potential will be larger and, correspondingly, the density of resonance 
curves and avoided crossings in Figs. E] and El will be increased too. This should increase 
the DR cross section (more states to be captured to) and the ratio Pd/ Pion (more avoided 
crossings) . 

Another approximation is that the n = 2 Rydberg states have dissociative character and 
are poorly described by this model. However, we believe that the dominant DR pathways 
are triggered by an initial capture into higher Rydberg states, rather than being directly 
captured into n = 2 states. 



VI. COMPARISON WITH OUR PREVIOUS THEORETICAL STUDY 



In the previous theoretical study of DR in HCO + published by some of us, 3l| the Renner- 
Teller coupling was not taken into account because it appeared to have a small effect on 
the potential surfaces of excited molecular states of the HCO molecule. However, Renner- 
Teller coupling involves degenerate ir states and therefore should play an important role in 
electron-ion scattering especially when pit and pa molecular potentials approach closely to 
each other. Although Renner- Teller coupling was not included in the previous study, some 
nonadiabatic effects were considered there. The quantum defects for sa,pa,da states 

obtained from ab initio calculation demonstrate sharp variations at some configurations Q 
of the three HCO nuclei. Therefore, instead of using the three adiabatic defects fJ,(Q), a 
numerical diabatization procedure was applied to give a Q-dependent 3x3 matrix ^ d \Q) 
of quantum defects. The Q dependence of the matrix is weaker and the matrix has more 
information about nonadiabatic effects than the diagonal matrix fi(Q). On the other hand 
the nonadiabatic effects represented by non-diagonal elements of /i(Q) might have a smaller 
effect on the dynamics of the electron comparing with the nonadiabatic couplings involving 
degenerate n states. Indeed, the DR cross section obtained in the previous study j^j] is 
much smaller than the cross section of the present treatment. To verify that the increase 
in the DR cross section is indeed due to the Renner- Teller coupling, we artificially set the 
non-diagonal elements in the Hamiltonian of Eq. to zero and calculated again the cross 
section for M — 0. The resulting DR cross section is about a factor of three smaller than the 
cross section shown in Fig. El and on average it is close to the result of the previous study. 
It is worth mentioning that the calculation of Ref. j^l] accounts for all three vibrational 
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degrees of freedom but the very high vibrational levels that generate the steep Rydberg 
n = 2 curves in Fig. |S] (which have large widths and therefore lead to an increase in the 
capture probability) were omitted; these two effects would tend to balance each other, and 
might give roughly the same cross section as in the present study when off-diagonal couplings 
are switched off. We speculate that if the CO bond vibration was included in the present 
treatment the effect of the Renner- Teller effect would likely be even more pronounced. To 
summarize our discussion of this point, we conclude that the Renner- Teller effect appears to 
be the most important nonadiabatic coupling in the dissociative recombination of HCO + . 

The theoretical cross section obtained in this study is inversely-proportional to the inci- 
dent energy below 0.1 eV (see Fig. EI). Above 0.1 eV the cross section falls more rapidly 
with increasing energy. This is caused by the fact that at 0.1 eV an additional ionic channel 
becomes open. That ionic channel is responsible for the large probability to capture the 
electron. Above that ionization threshold, the density of Rydberg states (in which capture 
can occur) in much lower, as can be seen from Fig. EJ and therefore the capture probability 
drops. Moreover, resonance states above that ionization threshold have another competing 
decay channel for autoionization, which steals flux away from the DR observable. For the 
ground vibrational level {OO ^} (dot-dashed line in Fig. |3]), that additional open channel 
is {Ol 1 ^}, for the first excited level there are two nearby levels: 102"^} and {02 2 t>3}. A 



similar behavior of the cross section was observed in DR of Hg 
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121: the cross section 



drops down once a new ionic channel becomes open. The cross section from the previous 
theoretical study of HCO + (the solid line in Fig. EJ) does not show that behavior. It behaves 
smoothly as the energy crosses the new ionization threshold. This is because the Renner- 
Teller effect was not accounted for: hence transitions of the type {OO ^} < — > {Ol 1 ^} are 
forbidden by the different symmetries of those vibrational states. 

Another difference with the previous study is the manner in which the final cross section 
is calculated. In Ref. [^, the scattering matrix for electron-ion collisions was explicitly 
calculated, but only in the electron-ion subset of channel space with no channel indices 
explicitly referring to dissociation. All of the channels \i) are vibronic states. However, 
some of the channels \i) are open for dissociation, and as a result, the Sa> matrix is not 
completely unitary. The defect from unitarity was used to calculate the dissociative flux and 
the corresponding DR cross section. The resulting cross section from this method has a rich 
structure with many Rydberg series of resonances. The cross section was then averaged over 
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an appropriate thermal distribution of electron energies. The averaged cross section still has 
a number of resonances that are not seen in the comparatively low resolution experiment. 
In the present paper we have employed a different averaging procedure, as was discussed 
above. This procedure gives a very smooth curve for the cross section. 

VII. CONCLUSION 

In summary, we would like to stress the following: 

• We have developed a theoretical approach to obtain the cross section for dissociative 
recombination of HCO + . The approach is based on a model reaction matrix K that in- 
cludes the non-Born-Oppenheimer Renner- Teller coupling between doubly-degenerate 
vibrational modes of HCO + and the degenerate continuum states of the incident elec- 
tron. The structure of the matrix K is derived directly from the Hamiltonian that 
accounts for the Renner- Teller physics, with elements of the constructed .fT-matrix, if 
diagonalized, that reproduce the quantum defects obtained from ab initio calculations. 

• We have shown that even a modest Renner- Teller splitting can cause a greatly en- 
hanced DR rate. The previous theory without the Renner- Teller coupling gave a lower 
DR cross section. The present DR cross sections are still generally about a factor of 2 
lower than the experimental cross sections. 

• The vibration along the CO bond was not included in the present model. We believe 
that the inclusion of the CO bond should increase the density of HCO states available 
to be populated during the electron-ion collision. This will probably increase the DR 
cross section, but the quantitative amount remains to be determined by future studies. 
Another error expected to be of order 70% is introduced by our use of an adiabatic 
approximation in the Rqh coordinate, and this should also be improved in future 
calculations. 

• In the present treatment the dissociation occurs predominantly through the indirect 
pathway, i.e. through Rydberg states of HCO. The estimated autoionization probabil- 
ity is comparable with the predissociation probability, but it seems likely that inclusion 
of CO vibrations may increase the probability of predissociation. 
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• The approach is quite general and can be applied to other closed-shell linear molecular 
ions, such as HCS + . 
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